m1=2500; m2=320; k1=80000; k2=500000; b1 = 350; b2 = 15020;
num1=[(m1+m2) b2 k2];
den1=[(m1*m2) (m1*(b1+b2))+(m2*b1) (m1*(k1+k2))+(m2*k1)+(b1*b2) (b1*k2)+(b2*k1) k1*k2];
'G(s)1'
printsys(num1,den1)
t=0:0.01:50;
y1= step(num1,den1,t);
plot(t,y1)
grid
title('open-loop response to unit step actuated force (transfer function)')